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As an unstable light pure leptonic system, positronium is a very specific probe atom to 
test bound state QED. In contrast to ordinary QED for free leptons, the bound state 
QED theory is not so well understood and bound state approaches deserve highly accu- 
rate tests. We present a brief overview of precision studies of positronium paying special 
attention to uncertainties of theory as well as comparison of theory and experiment. We 
also consider in detail advantages and disadvantages of positronium tests compared to 
other QED experiments. 
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1. Introduction 

Quantum electrodynamics (QED) is the only quantum field theory which can be 
successfully applied to a broad range of effects (bound states, scattering, decay) and 
energies from microwave radiation to high energies in the GeV range and deliver us 
various accurate predictions for measurable quantities with an uncertainty reaching 
the ppt level. 

QED of photons and leptons is an absolutely correct theory in a sense that its 
Lagrangian is well defined and in principle there is no problem for performing any 
calculations. 

• However, that is not sufficient since exact calculations are strongly limited by 
increasing difficulties in the calculation of higher-order effects and we always 
have to deal with a finite number of terms in the perturbative expansion. A 
question therefore arises how to estimate terms, which are too complicated to be 
calculated and sometimes that is in part art. 

• However, QED is in some way incomplete since electromagnetic interactions may 
involve hadrons and strong interactions which cannot be calculated ab initio. The 
weak interactions may also be involved, but usually it is not a problem to find 
related contributions. 

• However, there is a basic theoretical problem while performing comparison to 
experiment. Theory is not in position to predict any numbers. What theory can 
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only do is to express some measurable quantity in terms of others. In particular, 
to produce any quantitative prediction within bound state QED, we need to 
be first able to determine with a proper accuracy values of basic fundamental 
constants (the Rydberg constant Roc, the fine structure constant a as well as 
different masses and magnetic moments) and auxiliary parameters (such as the 
proton charge radius) due to the hadronic sector, which are necessary input data 
for bound state QED calculations. 

• However, there is no way to solve the bound state problem in general. The ap- 
plication of QED to the bound state forms a field called bound state QED, which 
experiences its own difficulties, additionally to QED problems for free particles. 

2. Bound state QED 

Bound state QED is quite attractive as a training field to solve the bound state prob- 
lem in quantum theory. It may be helpful for few-nucleon nuclei and for hadronic 
particles. In particular, there is a similarity in physics of positronium and quark- 
antiquark systems (mesons). 

The difficulties of precision QED calculations for free particles are mainly due to 
an increasing number (up to one thousand) of complicated diagrams (the four-loop 
level). The bound state QED theory deals with much simpler diagrams, however, the 
charged particles are bound there rather than free, and thus the Coulomb exchange 
may be not a small effect. A detailed review on QED calculations for light atoms 
can be found in^. 

The free QED involves only one small parameter a, while the bound state QED 
theory needs at least three and all three expansions are not quite goocH 

• Indeed, we still have to deal with a, the power of which indicates the number 
of QED loops involved. The expansion is asymptotical but that is not important 
since the bound state calculations mainly need one-loop and two-loop contribu- 
tions, with three-loop effects being important rather seldom. 

• The Coulomb strength Za appears because of binding effects and the parameters 
a and Za behave in a quite different way. There is a number of contributions 
where we need to sum over an infinite number of Coulomb exchanges (e.g. the 
Bethe logarithm) . Importance of all the exchanges assumes their essentially non- 
relativistic behaviour and still allows a Za expansion, which cannot be avoided 
since calculations exact in Za are possible for a few contributions only. However, 
that is not a well behaving expansion, because the limit Za — > is related 
to an unbound two-body system and thus leads to a non-analytic behaviour 
of perturbative expressions. The non-analyticity in Coulomb systems is usually 
accompanied with numerous logarithmic factors. For Z = 1 (hydrogen, muonium, 
positronium) one can find: ln(l/Za) ~ 5. The corrections known up to now may 
include up to cube of this logarithm, \n 3 (l/Za) ~ 120 (for the Lamb shift), and 
up to logarithm squared for hyperfine structure and positronium physics**!. For 
several reasons the non-logarithmic contributions also involve big coefficientJ^. 
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• A bound state problem supposes that we deal with an atom consisting of an 
orbiting particle (s) and an attractor (nucleus) and that involves one more pa- 
rameter, a ratio m/M of masses of the orbiting particle (mainly an electron) and 
the nucleus which for conventional atoms is m/M ~ 0.5 ■ 10~ 3 A _1 , for muonium 
is m/M ~ 1/207 and for positronium is m/M = 1. The behaviour of the expan- 
sion in m/M is not good since the limit m/M — > is related to a kind of bound 
"neutrino" . The non-analytic behaviour shows itself in logarithmic terms and, 
e.g., for muonium ln(M/m) ~ 5.5. 

• Some more parameters are involved due to nuclear effects, such as, e.g., a ratio 
of the Bohr radius to the nuclear radius. 

Thus, the bound state QED theory involves a rather rich spectrum of problems 
and it deserves to be tested, particularly in spite of the lack of well established uni- 
versal prescriptions appropriate for the two-body bound state problem in general. 
There are two basic problems of the precision bound state QED theory: Lamb shift 
and hyperfine structure. 

• The Lamb shift calculations mainly need an external field approximation, while 
recoil corrections are less important and only the simplest of them are involved. 

• In contrast, calculations of the hyperfme interval are crucially affected by the 
recoil effects, while some external field effects (e.g., the higher-order two-loop 
corrections) are relatively less important. 

3. Hyperfine structure in bound state QED 

Most of interest to positronium properties is due to recoil effects which are crucial 
since m/M = 1. For this reason we consider here in more detail studies of the 
hyperfine structure in light atoms. 

Magnetic effects are relativistic effects and thus, in contrast to the Coulomb 
interaction responsible for the Lamb shift, the higher momentum transfers and 
shorter distances are more important. At shorter distances the recoil and nuclear- 
structure effects arc enhanced. The nuclear effects in hydrogen and other light atoms 
dominate over the bound state QEbSG3 Still there are three possible QED tests 
with the hyperfine structure in which the problem of nuclear effects can be avoided. 

• A comparison of the Is and 2s hyperfine intervals that offers a specific differenceSE 

D 2 i = 8^hfs(2s)-^hfs(1s) , (1) 

which is immune to leading effects of the nuclear structure. 

• A comparison of conventional and muonic atoms for the same nucleus. 

• A study of a pure leptonic atomic system such as muonium and positronium. 

Determination of the 2s hyperfine interval in muonic hydrogen is in part the goal 
of an PSI experiment which is now in progress. The HFS interval in the Is and 
2s states was successfully studied for several light atoms. The more complicated 
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measurement is related to the metastable 2s stat cr^M Mi lz \ J) 21 theory is com- 
pared to experiment in Fig.^ The theory of D21 for several atoms is summarized in 
Tabled The dominant uncertainty of QED theory is due to higher-order one-loop 
and two-loop corrections in order a(Za) 3 and a 2 (Za) 2 and recoil corrections in 
order a(Za) 2 (m/M) (in units of the Is HFS splitting). The higher-order nuclear 
effects also substantially contribute to the uncertainty. 





optical, 2003 


rf, 2000 




rf, 1956 

Theory 



49.0 50.0 
Value of D 2 i in hydrogen [kHz] 




-1192.00 -1191.00 -1190.00 -1189.00 
Value of D 21 in 3 He + ion [kHz] 



Fig. 1. Determination of the D21 difference in hydrogen and helium-3 ion. The references can be 
found in Table 10. 



Table 1. Theory of the specific difference D21 in light 
hydrogen-like atoms 5,4 . 



Contribution 


Hydrogen 


Deuterium 


a He+ ion 




[kHz] 


[kHz] 


[kHz] 


D 2 i(QED3) 


48.937 


11.305 6 


-1189.252 


£>2i(QED4) 


0.018(3) 


0.004 3(5) 


-1.137(53) 


r>2i(micl) 


-0.002 


0.002 6(2) 


0.317(36) 


D 2 i(theo) 


48.953(3) 


11.312 5(5) 


-1190.083(63) 



The muonium theory of the Is hyperfine interval is summarized in Tabled The 
dominant uncertainty of QED theorjEH is due to higher-order recoil corrections 
in order a(Za) 2 (m/M) and (Za) 3 (m/M) which are also in part responsible for 
the uncertainty of D21 (see above). The other important part of the uncertainty is 
related to the determination of the leading term (so-called Fermi energy) 

^= 1 4a 2 ^R ocC (^-) 3 (2) 
h 3 hb \ rrifj. + m e J 



in terms of fundamental constants (i^/hb and a and inaccuracy in their determi- 
nation. 
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Table 2. Theory of the Is hy- 
pcrfinc splitting in muonium. The 
calculations 13 have been adjusted 
to a' 1 = 137.035 998 76(52) 14 and 
Mm/Mp = 3.183 34517(36) which was 
obtained from the analysis of the data on 
Breit-Rabi levels in muonium 15 and other 
less accurate experiments. 



Term 


AE 




[kHz] 


Ep 


4459 031.88(50) 


CLe 


5170.926(1) 


QED2 


- 873.147 


QED3 


- 26.410 


QED4 


- 0.551(218) 


Hadronic 


0.240(4) 


Weak 


- 0.065 


Total 


4 463 302.73(55) 



4. Why positronium ? 

As mentioned before, the recoil effects are better seen in positronium. Below we 
consider the positronium spectrum and a comparison of theory to experiment. We 
find that the uncertainty in calculating all experimentally studied transitions is 
related to the same recoil contributions as for the muonium HFS interval. 

The recoil effects play a crucial role in a two-body bound problem showing how 
closely the bound system is to a real two-body system. However, the significance of 
positronium is not limited by the possibility to verify the theory of recoil corrections. 

• Since the corrections of interest are enhanced (m/M is not a suppressing factor 
any longer), the fractional accuracy for successful high-precision tests is now 
relatively low. As a result, in contrast to hydrogen, the interpretation of the 
measurements of the IS — 2S interval does not crucially involve knowledge of 
the Rydberg constant with high accuracy. A study of the hyperfine interval does 
not require a value of the fine structure constant with high accuracy as it is 
in muonium. Since m/M = 1, it is not necessary to determine either m/M or 
A'Nuci/mb in an additional experiment. In other words, positronium offers several 
high precision tests of bound state QED without determinations of fundamental 
constants with high accuracy. 

• The hyperfine and recoil effects are enhanced and thus can be seen not only in a 
direct study of the hyperfine structure but also in the investigation of the gross 
or fine structure in contrast to hydrogen and muonium. Thus, positronium offers 
a few transitions which can be studied with high accuracy (see Table EJ). Adding 
to that an opportunity of different experiments on positronium annihilation, we 
find a big variety of properties to be studied. 

• It is not even necessary to mention that as a light pure leptonic atomic system, 
the positronium atom is free of hadronic effects. It is important that positronium 



1, 2008 22:13 WSPC/Guidelines-IJMPA preprint 



6 Savely G. Karshenboim 

is light because hadronic effects in leptonic systems involve a high momentum 
transfer. They can be seen in muonium and in the muon anomalous magnetic 
moment they are responsible for a dominant part of the theoretical uncertainty, 
while they are strongly suppressed for positronium and the anomalous magnetic 
moment of electron. 

Since positronium is a very specific atom, the notation for positronium is slightly 
different from other two-body atoms. First, in two-body atoms, even in hydrogen 
and muonium, it is customary to keep the value of the nuclear charge Z in order to 
recognize exchange photons and photons of QED radiative effects and thus to trace 
the origin of different corrections. In positronium there is an interference between 
both kinds of photons because of annihilation diagrams and thus it is meaningless 
to keep Z. 

Hydrogen and other two-body atoms are one-electron systems and one can use 
for them both small (e.g., Is) and capital (e.g., IS) letters to denote levels. The 
former are used for an electron, while the latter are for all electrons in an atom 
and that is indeed the same for single-electron atomic systems. Here, we prefer to 
use small-case letters for hydrogen and others. Since the nuclear spin effects are not 
suppressed, positronium has a structure of energy levels (in respect to their spin and 
angular momentum) rather similar to a two-electron system (such as the neutral 
helium atom) and one has to use only capital letter for its orbital momentum. 



Table 3. Theoretical predictions for positronium. 
Comparison to experimental data is presented in the 
figures quoted in the last column of the table. 



Quantity 


Prediction 


Figure 


Al/ H FS(1S) 


203 391.7(6) MHz 


m 


Au(l 3 Si -2 3 Si) 


1233 607 222.2(6) MHz 


m 


Ai/(2 3 Si - 2 3 P ) 


18 498.25(9) MHz 


□ 


Au(2 3 Si - 2 3 Pi) 


13 012.41(9) MHz 


m 


A;;(2 3 Si - 2 3 P 2 ) 


8 625.70(9) MHz 


m 


Ai/(2 3 Si - 2 1 P 1 ) 


11185.37(9) MHz 


m 


r(p-Ps) 


7989.62(4) ^s" 1 


m 


r(o-Ps) 


7.039 96(2) ims- 1 


c> 


Br 47/27 (pPs) 


1.439(2) ■ 10~ 6 


m 



We summarize theoretical predictions for various transitions and decay rates 
in Table 3. The Table does not contain experimental results but only references 
to figures where a comparison of theory to experiment is performed. Our results 
for decay rates and the HFS interval are slightly different from those quoted in 
the literature because we include into the decay rates contributions related to the 
decay to non-minimal number of photons (four for parapositronium and five for 
orthopositronium) which are sometimes omitted and because of our conservative 
estimation of uncertainty. To be conservative, we estimate the final uncertainty 
by half of the value of the leading logarithmic terms, if the non- leading term is 
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unknown or reduces the logarithmic contribution, and by half of the value of the 
whole logarithmic contribution, if the non-leading term enhances the leading term 
(cf. [16]). 

5. Positronium spectrum 

Let us now consider the positronium spectrum in more detail. The hyperfine inter- 
val in positronium has been measured with the highest absolute accuracy among 
positronium transitions. A comparison of theory and experiment is presented in 
Fig. [21 while theory is summarized in Table 0] The theoretical contributions are 
classified in units of the Fermi energy 

Ep = 12 c ^ m( ^ ' ( 3 ) 
which is defined including the virtual one-photon annihilation. E.g., QED2 is related 
to corrections of relative order a 2 . This notation is helpful for a comparison to the 
theory of decay rates, while the absolute units, also presented in the Table, allow a 
simple and direct comparison with theory of the IS* — 25 interval. Agreement with 
experiment is not perfect (within 2.5-3 standard deviations) which demonstrates 
that positronium is not very well understood. 

1 1 1 1 1 1 1 1 1 1 1 

a h 



b i * 

Theory 

i i i i i i i i i i i 

203 385 203 390 203 395 

Hyperfine interval in positronium [MHz] 

Fig. 2. The positronium IS hyperfine splitting: a comparison of theory to experiment (a - [17], 
6- [18]). 

The IS* — 25 interval was measured in orthopositronium by a method of three- 
photon ionization which has a resonance related to the two-photon excitation of 
the metastable 25 state from the ground state. The absolute uncertainty of the 
15 — 25 interval is compatible with the hyperfine splitting although somewhat 
lower. The uncertainty of the theoretical predictions for both hyperfine interval 
and the 15 — 25 transition (see Table [SJ) is determined by a correction of order 
a 7 mc 2 . That is related to a 3 Ep in the positronium hyperfine interval (QED3) 
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Table 4. QED contributions to the 15 hy- 
perfine interval in positronium. The uncer- 
tainty is presented following [16]. 

Term AE 
[MHz] 

E F oFmP 204 386.6 

QED1 a 5 mc 2 -1005.5 

QED2 a 6 mc 2 

QED3 aW . ^(ePZEBa 

Total 203 391.7(6) 



and thus to QEDA corrections in atoms with heavy nucleus (m/M <C 1), namely, 
a(Za) 2 (m/M)Ep and (Za) 3 (m/M)Ep. A comparison of theory to experiment is 
presented in Fig. [21 



Table 5. Theory of the 
l 3 5i — 2 3 Si interval in positron- 



ium. 



Term AE 

[MHz] 

~a I mc 1 1233 690 735.1 
a 4 mc 2 -82 005.6 

a 5 mc 2 -1 501A 

a 6 mc 2 -7.JSI 

a 7 mc 2 1.2(6P1 

Total 1233 607 222.2(6) 



a i — . — i 






Theory 





1 233 607 200 1 233 607 220 1 233 607 240 

Positronium 1S - 2S interval [MHz] 



Fig. 3. Theory and experiment (a - [25], b - [26]) for determination of the l 3 Si — 2 3 Si interval 
in positronium. 
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The theoretical uncertainty for the n = 2 fine structure in positronium (see 
Fig.Q is compatible with that for the hyperfme interval and the IS — 2S transition, 
but the experiment on 2S — 2P tran sitio ns cannot provide competitive results (see 
Fig. EJ. However, progress is possible^! 




orthopositronium (F=1) parapositronium (F=0) 



Fig. 4. Fine structure in positronium. The scale does not allow to see a width of all levels, but 
only of the broad 25 singlet state. The 2 3 S± — 2 1 Pi transition is forbidden, however, it may be 
observed by applying a magnetic field. 




8620 8630 8640 



11 170 11 180 11 190 13 000 13 010 13 020 

Positronium 2S-2P fine structure [MHz] 



18 500 18 510 18 520 



Fig. 5. Measurements of the fine structure transitions at n = 2. Experiments were performed at 
Brandeis University (a - [27]), University of Michigan (b — [28] and d — [30]) and Mainz University 
(c - [29] and e - [31]). 



6. Positronium annihilation 

Positronium is an unstable atom because of electron and positron annihilation. 
The IS state of positronium decays mainly into two or three photons depending 
on its spin. Thus, the orthopositronium (S = 1) annihilates into an odd number 
of photons, while the parapositronium (S = 0) into an even number. The leading 
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rates are 



lf>Ps) 



2(tt 2 -9) 6 toc 2 



9tt 



ri?( P Ps) 



, mc 



2 a — 



(4) 
(5) 



The lifetime of orthopositronium was responsible for a long-standing discrep- 
ancy of theory and most accurate experiments performed at the University of 
Michigar F^^ I while theory was in fair agreement with less accurate experiments, 
which disagreed with the Michigan data. The present situation is summarized in 
Fig. El Hopefully, the crisis seems to be resolved since the new Michigan vac- 
uum resulj2Hl now agrees with theory and the reevaluatiorP^I of the former gas 
experimental shifted the value towards theory. A further reexamination of system- 
atic effects is in progress and the final uncertainty will probably be higheJ^". 

Meanwhile, the lifetime of parapositronium also measured at the University of 
Michiga rP is in fair agreement with the theoretical prediction (see Fig.0). The the- 
oretical predictions for the annihilation decay rates for ortho- and parapositronium 
are summarized in Table El 




7.025 



7.035 7.045 
Decay rate of orthopositronium [u.sec-1 ] 



7.055 



Fig. 6. Determination of the orthopositronium annihilation decay rate: a comparison of theory 
to experiment (a - [34], b - [35], c - [32], d - [36], e - [33], / - [37], g - [38], h - [39]). 



We do not consider here exotic decay modes related to possible new physics, 
but rather rare modes possible within QED and the Standard Model. Two rare 
QED modes are allowed at a detectable level. Their theory and experiments are 
summarized in Table [7| We also include in the Table a photonless annihilation of 
orthopositronium into a pair of neutrino-antincutrino. The other modes with a pair 
of vv are also possible but their branching ratios are a few orders of magnitude 
lower. 
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Ann Arbor, 
vacuum 



-•—I 



Theory 



7980 7990 8000 

decay rate of parapositronium [u.sec _1 ] 





i ■ i 


I • 


— 1 e 


d l • 












I — 


• — I c 






h- 


b 




1 


a l • — 




Theory 
i.i.i 


1.0 


1.5 2.0 



Parapositronium branching 47/27 [ppm] 



Fig. 7. Annihilation of parapositronium: decay rate measured in University of Michigan in Ann- 
Arbor (left) and branching Br(47/27) measured around the world (a - [53], b - [54], c - [55], d - 
[56], e - [48]). 



Table 6. Theory of annihilation decay rate of ortho- and 
parapositroium (the IS state). The leading contributions are 
defined above in Eqs. (4) and (5). The decay rate of or- 
tho/parapositronium into five/four photons is included into 
QED2 terms. 



Contribution 


Decay rate 


Decay rate 




of orthopositronium 


of parapositronium 




[/isec ] 


[/jsec ] 


r(°) 


7.211 17 


8 032.50 


a ■ r(°) 


-0.172 30 


-47.25 


a 2 ■ r(°) 


0.001 ll(lpO 


4.43(ljSl 


Q 3 . r (o) 


-0.000 01(2p44| 


-0.08(4J3EI1 


Total 


7.039 96(2) 


7989.62(4) 



Table 7. Rare positronium decay modes. The theoretical uncertainty for 47 decay 
is related to unknown higher-order corrections. 



Rare mode 


Leading 


mode 


Branching (theor) 


Branching (exp) 


oPs - 


^5 7 


oPs - 


3 7 


1.0 -10- 6 , Ref.BSESl 


(2.2ff b s ) ■ IO" 6 , Ref.EZl 
1.7(11) ■ 10" 6 , Ref.EHl 


oPs - 


-» vv 


oPs - 


3 7 


6.2 -10- 21 , Ref.EHl 


< 5.8 • 10" 4 (90% CLPJ" 

< 2.8 • 10" 4 (90% CLpU 


pPs - 


-> 47 


pPs - 


.27 


1.439(2) • 10" 6 Ref.l^l 


1.30(30) ■ 10" 6 , Ref.ElH 
1.50(11) ■ 10" 6 , Ref.ES 
1.48(17) ■ 10" 6 , Ref.EHl 
1.19(26) ■ 10" 6 , Ref.EHl 
1.14(39) ■ 10" 6 , Ref.EHl 



We have to note that the branching fraction of the orthopositronium decay into 
neutrinos is so low, that this decay is unlikely to be soon detected, however, a mode 
oPs — * nothing (since the neutrino is not detectable) can still be of interest if the 
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neutrino has a non- van ishi ng magnetic moment (see, e.g., Fig. |SJ). This problem 
was discussed in part in ISZ but we do not like to present here any numbers. In our 
opinion, results of such analysis can depend on a model introducing neutrino mass 
and magnetic moment. Still, studies of the pure neutrino modes can provide a limit 
on the r-neutrino magnetic moment at a level above (but possibly not much above) 
the current limits^! Since systematic effects are different and the interpretation 
depends on the model, we think it may be important to have several independent 
limitations. 



V V 




V 



Fig. 8. Annihilation of orthopositronium into a pair of uv: via Z-boson (left) and photon (right). 
The former diagram is related to the Standard Model, while the latter is present only if /i„ ^ 0. 



7. Positronium and other QED tests 

Tests with positronium are quite different from other tests. Advantages and dis- 
agvantages of positronium studies and some other QED experiments are listed in 
Table El 



Table 8. Advantages and disadvantages of different QED tests for conventional light atoms, 
pure leptonic bound systems and free leptons. 



System 


Value 


Dominant uncertainties 


Hydrogen 


Lamb shift 


Nuclear size, higher-order two-loop effects, Roo 


Hydrogen 


Is HFS 


Huge uncertainty due to nuclear structure 


Hydrogen 




Experiment 


3 Hc+ 


£>21 


Higher-order two-loop, recoil, nuclear effects, experiment 


4 He+ 


Lamb shift 


Nuclear size, higher-order two-loop effects, experiment 


Muonium 


Is HFS 


Higher-order recoil effects, Hb, oc 


Positronium 


15 HFS 


Higher-order recoil effects, experiment 


Positronium 


IS - 2S 


Higher-order recoil effects, experiment 


Positronium 


2S - 2P 


Experiment 


Orthopositronium 


r(la) 


Experiment 


Parapositronium 


r(is) 


Experiment 


Electron 


9-2 


Uncertainty: a, cavity QED effects 


Muon 


9-2 


Uncertainty: hadronic effects 



To compare different QED tests with/without positronium, we also check which 
contributions are crucial for a comparison of theory to experiments. They are sum- 
marized in Table |5J 
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Table 9. Comparison of the bound state QED theory and experiment: crucial orders of magnitude 
for the energy levels and decay rates needed for the comparison 2 (in units of mc 2 ). 

Atom, value Crucial order(s) 

[mc 2 ] 

Hydrogen (gross structure) a(Za)' , OL 1 (Zaf' 

Hydrogen (fine structure) ct(Za) 7 , a 2 (Zct) f ' 

Hydrogen (Lamb shift) ', a 2 (Za) 6 

3 He+ ion (2s HFS) a(Za) 7 m/M, a 2 (Za) 6 m/A4, a(Zaf(m/M) 2 , (Za) 7 (m/M) 2 

4 Hc+ ion (Lamb shift) a(Za) 7 , a 2 (Za) 6 

Muonium (Is HFS) a(Za) 7 m/M, a(Z a) 6 (m / M) 2 , (Za) 7 {m/M) 2 



Positronium (IS HFS) 


a 1 


Positronium (IS - 23) 


a 7 


Positronium (2S - 2P) 


a 7 


Para-positronium (decay rate) 


a 7 


Ortho-positronium (decay rate) 


a 8 


Para-positronium (47 branching) 


a 8 


Ortho-positronium (57 branching) 


a 8 



8. Higher-order logarithmic corrections and uncertainty of 
positronium calculations 

To consider HFS tests in detail, we summarize the results on a study of the hyperfme 
structure in light two-body atoms in Table ITU1 In these tests the QED uncertainty 
has never been a limiting factor for comparison of theory and experiment. The 
following uncertainties may also be involved: due to nuclear effects, determination 
of fundamental constants or measurement. One of them is always bigger than the 
theoretical uncertainty. 

The theory and experiment are in general in good agreement. The uncertainty 
for the Is hyperfme interval in muonium and positronium and in part of D21 for 
helium-3 ion is related to the same recoil corrections in order a 3 (m/Af (see 
Table EJ. The positronium uncertainty in fractional units is much higher than that 
for experiments with heavy atoms (m/M < 1), however, its sensitivity to higher- 
order recoil effects is at approximately the same level. 

The higher-order recoil contributions are known only in the logarithmic approx- 
imation and we estimate the uncertainty as half the value of the leading term if the 
next-to-leading term is unknown or cancels a part of the leading term. However, 
if the leading and the next-to-the leading terms are of the same sign, we estimate 
the uncertainty as a half value of the whole logarithmic contribution (cf. ESEEEJ). 
Thus, a calculation of the next-to-leading terms does not reduce an uncertainty. A 
reason for such a conservative estimation is that the leading logarithmic term is 
mostly a result of a single contribution without any cancellation. For the ns state 
it is state-independent. For these reasons it has a kind of natural value and can be 
used to estimate the non-leading terms. In contrast, the non-leading term has quite 
accidental value and may be sometimes quite below the natural level. Thus, it can- 
not be used alone to estimate properly the next term of the logarithmic expansion. 
However, we cannot simply ignore the next-to-leading term, particularly in the case 
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when it is not too small. In the present paper we follow a compromise which allows 
to achieve a conservative estimation of the uncertainty. A calculation of the next 
of the non-leading terms helps us to check the reliability of this estimation using 
only the leading term and meantime is a necessary step towards a calculation of the 
whole contribution beyond the logarithmic approximation which is strongly needed 
for muonium and positronium HFS. 

Table 10. Comparison of experiment and theory for hyperfine structure in hydro- 
gen-like atoms. In the D21 case the reference is given only for the 2s hyperfine interval. 



Atom Exp. Theory A/a cr/E F 

[kHz] [kHz] [ppm] 

Hydrogen, D 2 i 49.13(13/'" 48.953(3) LI 0.09 

Hydrogen, D21 48.53(23^1 -1.8 0.16 

Hydrogen, D 2 i 49.13(40^ 0.4 0.28 

Deuterium, D 2 i 11.16(16^1 11.3125(5) -1.0 0.49 

3 Hc+ ion, D21 -1 189.979(7l]ni -1190.083(63) 1.10 0.01 

3 He+, D 2 i -1 190.1(16pH O0 0.18 

Muonium, Is 4 463 302.78(5^ 4 463 302.88(55) -0.18 0.11 
Positronium, IS 203 389 100(740pU 203 391700(600) -2.9 4.4 
Positronium, 15 203 397 500(1600) 19 -Z5 8.2 



9. Soft and hard QED effects 

Not only essential orders of QED corrections offer a possibility of comparing the 
efficiency of different experiments. The bound state QED theory clearly recognizes 
two kinds of contributions: the soft-photon contribution and the hard-photon con- 
tribution. The latter are very similar to free QED, while the former essentially 
involve binding effects. There are two most important soft-photon contributions. 

• Crucial corrections in the external field approximation are due to the higher-order 
two-loop self-energy. The inaccuracy in its calculation determines an uncertainty 
of the Lamb shift calculations for hydrogen and hydrogen-like ions, while similar 
effects involving the magnetic field significantly contribute to the uncertainty of 
D21 in the helium-3 ion. 

• Recoil corrections are crucially important for tests with hyperfine structure and 
they determine the uncertainty of muonium HFS, positronium spectrum (hyper- 
fine interval, Is — 2s transition, fine structure) and a part of the uncertainty of 
D21 for the helium-3 ion. 

Effects of the hard-photon exchange are very similar to effects of free QED. We 
note, however, that the regions of integration in momentum space are different. The 
most accurate free QED calculations which may be compared to experiment are 
related to the anomalous magnetic moments of electron and muon. The integration 
for them is performed over a kind of isotropic region in Euclidean space (|fco| ~ |k|) 
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and the crucial level is the four-loop approximation. For the bound state problems 
there are two other specific regions of integration. 

• For some problems it is sufficient to apply an external field approximation to 
hard-photon corrections and thus loops of exchange photons are related to zero 
energy transfer (fco = 0). Even for recoil effects the integration is essentially 
not covariant in Euclidean space including a contribution from a specific region 
|feo| -C |k|. The highest crucial orders are related to four- loop corrections for 
the external-field approximation and three-loop corrections for recoil effects. In 
contrast to the anomalous magnetic moment, that is a calculation for two different 
particles with its own simplifications and difficulties. 

• The other specific situation for integration is related to the positronium anni- 
hilation when some photons are real (ko — |k|). The orthopositronium studies 
(main decay mode and branching fraction for five-photon decay) allow to check 
calculations of four-loop corrections in such a non-isotropic region of integration. 
The accuracy of the branching fraction determination for the four-photon decay 
of parapositronium is at the level of a few percent. Because of a relatively large 
a-contribution to the 47 annihilation rate"*' 

r(pPs -> 4 7 ) = r (0) (pPs -> 4 7 ) x (l - 14.5(6)-) (6) 

the experiments are approaching a level, where four-loop diagrams are important. 

Thus, we see that the bound state problems supply us with a possibility to check 
modern algorithms for four-loop calculations and that is competitive to the anoma- 
lous magnetic moment of electron. We note, however, that the difficulties which 
appear at the four-loop level are a large number of digrams with high-dimension 
integrations and overlapping UV divergencies. Neither annihilation nor exchange 
loops are ultraviolet divergent, but the exchange loops involve strong infrared di- 
vergency. 

One can also compare^! soft-QED effects of the virtual one-photon annihilation 
and the real two- and three-photon annihilation. Soft Coulomb corrections to the 
hard annihilation block can be presented in fractional units. A crucial theoretical 
uncertainty related to a 3 corrections (QED3 term) is at the level of few ppm, which 
can be hardly achieved for a measurement of the annihilation rates. 

10. Summary 

Thus, we summarize our paper with the following statements. 

• Positronium spectroscopy offers a reliable test of our understanding of higher- 
order recoil corrections within bound state QED, which play a significant role 
in other QED tests with the hyperfine structure. The theoretical uncertainty for 
positronium is of a pure QED origin and other effects such as the nuclear structure 
effects or inaccuracy in determination of fundamental constants are not involved. 
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Studies of recoil effects are of particular interest because their significance answers 
a question whether we study a really two-body system. In conventional atoms 
(hydrogen, helium) the role of the nucleus as a particle is reduced. The positron- 
ium is a truly two-body system which is the closest to the neutral helium atom. In 
contrast to two-body atoms, the QED uncertainty in helium is bigger than that 
of experiment and theory needs a significant improvement. Thus, positronium 
theory serves as an intermediate step between hydrogen and helium. 

• A study of positronium annihilation including rare decay modes allows to test 
approaches to the calculation of four-loop diagrams but for a specific integra- 
tion region. Today the four-loop approximation is the highest level for precision 
calculations of any measurable QED quantities (if indeed no big logarithms are 
involved such as those in scattering kinematics). 

• There are exotic decay modes of positronium beyond QED, but within the Stan- 
dard Model and in particular a decay of orthopositronium into a pair of neutrino 
and antineutrino, which is a dominant mode beyond QED. The branching frac- 
tion is very low and it cannot be detected presently. Still, a study of this channel 
can probably give a limit for the magnetic moment of the r neutrino which, al- 
though somewhat weaker than the current limits, will, however, have completely 
different systematic effects. 
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